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Abstract 

The stability is analyzed of the equilibrium configurations of a mag- 
netic junction with a free layer that has cubic symmetry and two anisotropy 
j^ . axes in the layer plane. Different variants of the switching between vari- 

ed ' ous configurations are considered. A possibility is shown of the substan- 

I , tial lowering of the threshold current density needed for the switching. 

'■^ ' Numerical simulation is made of the switching dynamics for various con- 

(Ih ' figurations. 

o ■ 

1 Introduction 

►^ . Switching magnetic junctions by a spin-polarized current is one of the 

r~^ ' main spintronic effects. Besides the academic interest, this phenomenon 

^^ , may be used for high-density information processing, since the character- 

("^ ■ istic scales are the exchange interaction and the spin diffusion lengths of 

O^ ' the order of tens of nanometers. 

Nowadays interest has revived to magnetic junctions the layers of 
(^ • which have cubic magnetic anisotropy; the well-studied iron may be an 

f^ ' example [Tl[51[3]. Thin Fe(OOl) films have two equivalent anisotropy axes, 

[100] and [010], in the layer plane. This allows switching the layer mag- 
netizations between different easy axes by means of magnetic field and/or 
spin-polarized current, which may be used in memory cells with more than 
two stable states. 
$_( ' In the present work, we consider a magnetic junction with cubic anisotropy 

5^ , of the free layer placed between pinned and nonmagnetic ones. Switching 

such systems by applied magnetic field has been studied in Refs. [2l|4]. 
Here we investigate the switching by spin-polarized current at various rel- 
ative orientations of the pinned and free layer magnetization vectors. The 
thickness of the free layer is assumed to be small compared to the spin 
diffusion length, so that the macrospin approximation is valid [5]. 
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Two main mechanisms are known of the interaction between spin- 
polarized current and magnetic lattice, namely, spin transfer torque (STT) [S] 
[7] and spin injection leading to appearing regions of nonequilibrium spin 
polarization near the interfaces [21 12] (the term "field-like torque" is used 
sometimes in literature, because the mechanism action is equivalent to 
influence of some effective magnetic field in some cases). The first of the 
mechanisms mentioned (STT) is related with appearing a negative damp- 
ing that prevails over the positive Gilbert damping above the threshold 
current density; this leads to instability of the original magnetic configu- 
ration. The other ("injection") mechanism is based on increasing the sd 
exchange interaction energy between the nonequilibrium conduction elec- 
trons and the lattice under spin injection to the free layer of the magnetic 
junction, so that the original state becomes unstable, and a reorientation 
phase transition occurs. A unified theory of switching magnetic junctions 
including actions of both mechanisms was presented in Refs. [10U11| . The 
relative contribution of the mechanisms indicated depends on the layer 
parameters and the applied magnetic field. Earlier, the conditions were 
formulated |12II131 [5] under which the injection mechanism plays the main 
role (see below for details). In this work, special attention is paid to the 
latter. 

The main equations describing the magnetization of the thin free layer 
in the magnetic junction under spin-polarized current are presented in 
Sec. (21 In Sec. [3] we analyze the stability of the stationary configurations 
depending on the current density through the magnetic junction, as well 
as possible switching between these configurations. In Sec. |4] the current- 
voltage characteristic is found of the magnetic tunnel junction under for- 
ward and backward currents depending on the original configurations. In 
Sec. [5] the junction switching dynamics is simulated numerically. 

2 The main equations 

Let us consider a magnetic junction consisting of a pinned ferromagnetic 
layer 1, a free ferromagnetic layer 2, and a nonmagnetic layer 3 which 
closes the electric circuit. There is a thin spacer between the layers 1 and 
2 that prevents the direct exchange interaction between the magnetic lat- 
tices of the layers. The electric current flows perpendicular to the layers 
(CPP mode). The free layer 2 has cubic symmetry with three mutually or- 
thogonal symmetry axes and, correspondingly, magnetic anisotropy axes, 
one of which, [001], is perpendicular to the layer plane, while two other, 
[100] and [010], lie in the plane. The anisotropy energy of that layer (per 
area unit) is [14| 

Ua = ^MHaL\^ (m ■ ni)'' (m • na) V (m ■ na)"^ (m ■ ng) ' 

-H(M-n3)'(M-ni)'}, (1) 

where L is the thickness of layer 2, M is the saturation magnetization of 
that layer, M = M/|M| is the unit vector along the magnetization, ni, 
n2, ns are the unit vectors along [100], [010] and [001] axes, respectively. 
Ha is the effective anisotropy field. 



The free layer thickness L is assumed to be small compared to the 
spin diffusion length / and the inhomogeneity scale of the magnetic lattice 
in that layer (such a scale, the measure of the "spatial inertia" of the 
lattice, is the domain wall thickness 5). Under such conditions, layer 2 
manifests itself as united whole ( "macrospin" ) in respect of its magnetic 
behavior. This leads to the modification of the Landau-Lifshitz-Gilbert 
equation for the layer magnetization with the disappearance of the spatial 
derivative term and the introduction of a new term describing the current 
effects. With the cubic anisotropy taken into account, the equation takes 
the form (cf. 5 ) 
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Here Mi is the unit vector along the magnetization of layer 1, H is the 
applied magnetic field, Hd is the demagnetization field, k is the Gilbert 
damping factor, 7 is the gyromagnetic ratio, a = jHaS^ is the lattice 
magnetization diffusion constant. 
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where e is the electron charge, /is is the Bohr magneton, a is the sd 
exchange interaction constant, r is the spin relaxation time, A = L/l ^ 1, 
Q is the conduction spin polarization. 
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is the layer spin resistance [13], p is the electric resistivity; the quanti- 
ties without index refer to the free layer 2. The b = (aiMin) / (aMr) 
parameter describes influence of the pinned layer 1. 



The parameters p and k in the right-hand side of Eq. ([2]) describe in- 
jection and STT mechanisms of the spin-polarized current effect on the 
magnetic lattice under current flowing in the "forward" direction, corre- 
sponding to the electron drift in the f — > 2 — >■ 3 direction. A substitution 

(Ml ■ M) — > (Ml ■ M) , j ^^ ^|j| corresponds to the backward direc- 
tion (3 -^ 2 -i> f). 

As analysis shows [12], the intensity of the spin injection to the free 
layer is determined with the relationship between the spin resistances of 

the layers. At -^, -^ <^\ <^\ relationship between the spin resistances 

Z\ Zz 
of the pinned layer 1, free layer 2 and nonmagnetic layer 3, high spin 

injection level is reached due to effective injection from layer 1 to layer 2 

and "locking up" the injection from layer 2 to layer 3 13 . This leads to 

lowering the switching threshold due to the injection mechanism. Another 

possibility of the threshold lowering consists in applying external magnetic 

field [5] . The STT contribution becomes unimportant under low threshold 

of the injection switching magnetic junction, so that we will consider only 

the injection mechanism below. 

Under assumptions mentioned, the Landau-Lifshitz equation describ- 
ing the lattice dynamics of the free layer in the absence of damping has 
an integral of motion even in presence of a current through the junc- 
tion [151 I16| . Such an integral is the magnetic energy, including the 
Zeeman energy in applied magnetic field, the anisotropy energy, the de- 
magnetization energy, and the sd exchange interaction energy between 
conduction electrons and magnetic lattice. 

Let us consider a configuration with x axis along the current, yz plane 
parallel to the layer planes H = {0, .ffsin-i/i, Hcosip}, ni — {0, 0, 1}, 
na = {0, 1, 0}, ng = {0, 0, 1}, Hd = -47rM{M,, 0, 0}, Mi = {0, 0,1}. 

In spherical coordinates with the polar axis along [100] axis, M = 
{smO cosij)-, soi-OsirKJ), cos 6'}, the dimensionless (in MHaL units) mag- 
netic energy under forward current takes the form 

U(9, (j>) H /\,l-2/, 2/),l-4/,-2, 2, 

— ■ COS^t* — V) + o ^^^ '^^^ "*" o ^^"^ ^^^ 'r '"^^ 'r 
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The corresponding formula for backward current is obtained with sub- 
stitution 008 6^ — )■ (cos 6)"^, j — >■ — IjI in the last term of Eq. ([6]) describing 
the current effect: 

U{9, (j>) H . /\, 1-2/1 2a,l-4/i-2, 2, 

' = — — — cos(6' ~ y-') + - sm o cos 9 + - sm o sm cp cos (f> 
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Figure 1: The magnetic junction energy as a function of the angle between the 
magnetization vectors of the pinned and free layers with various values of the 
(dimensionless) forward current density. 



3 Stationary states and variants of switch- 
ing 

The stationary states of the system in study correspond to the extrema 
of the U{9, (j}) function; the minima of the function correspond to the 
stable equilibrium states. Because of the positive definiteness of the term 
with the azimuthal angle </>, it is sufficient to consider only the energy 
dependence on the polar angle 9 at fixed value (f> = 90° during the minima 
finding (this corresponds to in- plane position of the magnetization vector). 

There are three minima in absence of magnetic field {H = 0) and 
current (j = 0): 6 = 0° , 9 = 90°, and 9 — f80°, corresponding to the 
parallel, perpendicular and antiparallel relative orientations of the pinned 
and free layers. 

Let us present the stability analysis results for these stationary states 

in presence of the current under intense injection conditions -=-, -=- -C 

A ^ 1 (as mentioned above, the latter condition allows neglecting the 
STT contribution). The results are illustrated with plots of the U energy 
as a function of the 6 angle (Figs. [T]and[2]) at various values of the current 
density in forward and backward directions. 

1) The parallel configuration {9 — ) in absence of magnetic field 
{H = 0) is stable under forward and backward currents. 

2) Under the forward current increasing from zero, the energy mini- 
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Figure 2: The magnetic junction energy as a function of the angle between the 
magnetization vectors of the pinned and free layers with various values of the 
(dimensionless) backward current density. 



mum corresponding to the original perpendicular configuration {0 = 90°) 
shifts in the direction of the parallel configuration (i.e., 8 angle decreases). 
This takes place up to the current density value j = jo\/2/27 ~ 0.272 jo, 
when the deviation from the original position reaches arcsin (l/\/6) ~ 24° 
(i.e., 9 ~ 66° ). At this value, the energy minimum disappears (it changes 
to an inflection point), and the system switches abruptly to a parallel 
configuration and remains in that configuration under further variations 
of the current density, so that the 90° — >■ 0° switching is irreversible. 

Under backward current increasing from zero, the minimum corre- 
sponding to the original perpendicular configuration {6 = 90°) shifts in 
the direction of the antiparallel configuration, and the corresponding 9 
angle increases up to some value 6*0 (j) that tends to 
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at the high current limit. Under returning to zero current, the perpendic- 
ular orientation restores, so that the switching by the backward current 
is of "temporary" character. 

3) The antiparallel configuration {9 — 180°) in a magnetic field parallel 
to the magnetization of the pinned layer {tp = 0°) becomes unstable and 
switched to parallel one under high enough forward current. The threshold 



current density is 

^^•^ -'' 1 + iZ,/ZsX) iZ,/Z,X) V-H-J- ^» V^H-J- (') 

We see from Eq. © the mentioned possibility of the switching thresh- 
old lowering by means of an applied magnetic field close to (but lower 
than) the anisotropy field. Such an assistance of the magnetic field does 
not break the local character of the switching, because the magnetic field 
lower than the anisotropy field cannot do switching alone (without a cur- 
rent) . 

The parallel configuration appeared after switching is stable against 
further variations of the current, so that the switching by the forward 
current is irreversible. 

Under the backward current, the antiparallel configuration becomes 
unstable at the same (in magnitude) current density \j\ — jth, however, 
in this case switching takes place to a nonequilibrium stationary state 
9 = 6o(jth) or to the symmetrical state 9 — 360° — &o(jth)- With returning 
to zero current, the system does not return to antiparallel configuration, 
but comes to one of two perpendicular configurations 9 = 90° or 9 = 270° . 

Thus the following variants are possible of the switching between sta- 
tionary states: 1) switching an antiparallel configuration to a parallel one 
by turning up the forward current of j > jth density and subsequent turn- 
ing off; 2) switching an antiparallel configuration to a perpendicular one 
by turning up the backward current of the same density and subsequent 
turning off; 3) switching a perpendicular configuration to a parallel one by 
turning up and subsequent turning off the forward current of substantially 
lower density j > 0. 272 jo- 

4 Resistance of the magnetic tunnel junc- 
tion 

In experiments, the current-driven switching magnetic junction manifests 
itself, in the first place, as a change of the junction resistance. The re- 
sistance depends substantially on the relative orientation of the layers 
forming the junction; this is the cause of the well-known tunnel magne- 
toresistance effect. 

The conductance of a magnetic tunnel junction with 9 angle between 
the magnetization vectors of the layers takes the form |17) 

G(6») = Gp cos'' ^ + Gap sin' ^, (10) 

where Gp, Gap are the junction conductances at parallel {9 — 0°) and 
antiparallel (6^ = 180° ) relative orientation of the layers, respectively. 

It is convenient to describe the change of the junction resistance with 
the following ratio: 

R{9)~Rp _ p{l-cos9) 

Rp 2 + p{l + cos9)' ^ ' 
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Figure 3: The magnetic tunnel junction resistance as a function of the (di- 
mensionless) current density. The dots on the ordinate show the stationary 
configurations without current: parahel (square), perpendicular (circle), and 
antiparallel (skew cross). The arrows with corresponding tails show the resis- 
tance changes under the current change for different initial configurations. 
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Figure 4: Switching dynamics of the antiparallel configuration to parallel one 
by the forward current. 



where R{e) = 1/G{e), Rp = 1/Gp; p = [7?(180°) - Rp]/Rp is the tunnel 
magnetoresistance defined by usual way |18) . 

To find the resistance dependence on the current direction and density 
R{j), it is necessary to substitute 9{j) dependence to Eq. (|ll|l . With the 
foregoing analysis taking into account, the results are obtained shown in 
Fig. O A possibility is seen of the switching between different stationary 
states corresponding to different electric resistances. 



5 S imulat ion of t he magnet ic j unct ion s wit ch- 
ing dynamics 

Together with the investigation of the stationary states and the switching 
between them, the switching dynamics is of great interest, because it 
determines the speed of response of the devices based on the magnetic 
junctions. 

The time-dependent vector equation ((Sjl describing the dynamics with 
using polar coordinates {9, (j)) takes the form of a set of equations 
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Figure 5: Switching dynamics of the perpendicular configuration to parallel one 
by the forward current. 



where 
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The P{9) and K{9) functions describing the current effect take the 
form 
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Figure 6: Switching dynamics of the antiparallel configuration to perpendicular 
one by the backward current. The step shows the current turning-off time. 
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Figure 7: The perpendicular configuration evolution under turning on and turn- 
ing off a rectangular pulse of the backward current. 
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(19) 
for the bacfcward current. 

Since /la ^ 1 in fact, the precession of the free layer magnetization 
vector takes place in the layer plane on the whole, i.e., the vector end 
describes an ellipse elongated strongly in the plane. So the small terms 
proportional to ha cos^ (f> were omitted in derivation of Eqs. ((T2))-((T3}. 

In Ref. [16] an analytical solution was found of Eqs. (|12p . (|13|) at 
K = Q, h <gil, ha <Si ^ 'm the limiting cases of the zero damping (k = 0) 
and strong damping (k 3> 1). In the general case, it is reasonable to use 
numerical simulation. We used the Simulink software of the MATLAB 
system that is intended for simulating dynamical systems [19| . 

The following parameter values were given: k = 0.03, ha = 0.01, h = 
0, A = 0.1, Z2/ZX = Z2/ZS = 0.01, ar-^M = 60 (cf. [5]). The thermal 
noises initiating deviation of the free layer magnetization from the original 
unstable equilibrium state were imitated with giving a small initial devia- 
tion from such a state by an angle of 0.01 radians in the layer plane where 
the demagnetization field does not prevent fluctuation-induced deviations, 
so that minimal fluctuation energy is needed. The time dependences of 
the magnetization deviation from [100] axis 6{T) under forward current of 
j = 2jo density for antiparallel and perpendicular original configuration. 
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respectively, are shown in Figs. [4] and (5] the similar ones for the backward 
current are shown in Figs. [5] and [T] The dimensionless time is laid off as 
abscissa with to = (1 + k^)/(Atx^M) as the time unit; at M = 900 G one 
nanosecond corresponds to 200 scale divisions of the abscissa. The steps 
in Figs. [S] and [7] show the current turning on and turning off times. 

The numerical solution (simulation) results consist completely with 
foregoing analysis based on the angular dependence of the magnetic en- 
ergy. It is seen that direct switching occurs of the antiparallel and perpen- 
dicular configurations to parallel one under the forward current. Under 
the backward current, the switching occurs to an intermediate nonequilib- 
rium stationary state, from which a transition (or return, when the initial 
configuration is perpendicular) occurs to the perpendicular state. As was 
to be expected, the switching is accompanied with damped oscillation due 
to precession of the magnetization vector. 

With given parameter values, the characteristic switching times are 
of the order of nanoseconds, while the oscillation period is of the order 
of fractions of nanosecond. With increasing the magnetization and the 
damping constant, the speed of response rises (the latter up to some limits, 
because the switching process becomes aperiodic and slows with increasing 
damping at too strong damping (k ^ 1)). 

6 Conclusion 

The analysis shows a possibility of increasing the number of the switch- 
able states by using magnetic junctions with cubic-anisotropy layers. The 
fact is of interest that the switching of the perpendicular configuration 
to the parallel one requires current density lower by several times, than 
the switching of the antiparallel configuration. The fruitfulness should be 
noted of the combination of the energy approach to determining stationary 
states with numerical simulation of the switching processes. 
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